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Abstract The power spectrum is widely used in astronomy, to analyze temporal 
or spatial structure. In cosmology, it is used to quantify large-scale structure (LSS) 
and the cosmic microwave background (CMB). This is because the power spec- 
trum completely quantifies Gaussian random fields, which the CMB and LSS fields 
seem to be at early epochs. However, at late epochs and small scales, cosmological 
density fields become highly non-Gaussian. The power spectrum loses power to de- 
scribe LSS and CMB fields on small scales, most obviously through high covariance 
in the power spectrum as a function of scale. Practically, this significantly degrades 
constraints that observations can place on cosmological parameters. However, if 
a nonlinear transformation that produces a (more) Gaussian 1 -point distribution is 
applied to a field, the extra covariance in the field's power spectrum can be dra- 
matically reduced. In the case of the roughly lognormal low-redshift matter density 
field, a log transform accomplishes this. Applying a log transform to the density 
field before measuring the power spectrum also tightens cosmological parameter 
constraints by a factor of several. 

A Gaussian random field has convenient statistical properties. Its meaningful 
information is fully quantified by the power spectrum; all connected higher-order 
statistics vanish. Of particular importance for a measurer of (parameters which de- 
pend on) the power spectrum, all off-diagonal power-spectrum covariance matrix 
elements vanish for a Gaussian random field. 

On small scales at late epochs, the cosmological overdensity field 8 — — 1 is 
highly non-Gaussian, with high off-diagonal power-spectrum covariance. However, 
this non-Gaussianity is not of an essential form, and can be largely removed with 
a monotonic transformation. The near-lognormality (H of the density field may be 
exploited; a Gaussianizing transform like 8 — > ln(l + 8) much reduces the power- 
spectrum covariance [2|. As shown in Fig. [2} it also reduces error bars on inferred 
parameters, reaching a factor of 5 reduction in the best case of the tilt n s . 

Mark C. Neyrinck 

The Johns Hopkins University, Baltimore, MD 21218 



1 



2 



Mark C. Neyrinck 





o.l 

k [fc/Mpc] 

Fig. 1 Left: 8 in a 2-/i -1 Mpc slice of the 500-/i -1 Mpc Millennium simulation (MS), viewed 
with an unfortunate linear color scale. Middle: the same slice with a logarithmic color scale. Right: 
the 2D power spectra Pg and P\ n (\ + g) of S (black) and ln(l + S) (green), in 9 such slices. The wild, 
coherent fluctuations in Pg illustrate its high (co)variance, absent in Pi n (i+S)- 



10" 1 



*-< 1 

wicr 1 



10 



10 



.02 



ln(l +.5) 



^^^^ In al 















0.1 



0.1 







.05 



0.5 



Fig. 2 Fisher-matrix estimates of error-bar (half-)widths and error ellipses for the cosmological 
parameters lno"| and n s , inferred analyzing Pg (black) and fi n ( i+<5) (green). We show how they 
depend on the maximum wavenumber £ max included in a power- spectrum analysis of a 1-Gpc 
real-space matter density field. Diagonal panels show unmarginalized error bars over single pa- 
rameters. In the lower-left panel, errors in each parameter are marginalized over the other. The 
upper-right panel shows how error ellipses contract as £ nlax increases. There is an ellipse shown 
for each £ max constituting the curves in the other panels. Outside the bold ellipses, analyzing only 
large scales where £ max < 0.1 Mpc -1 , Pg and Pi n (i + a) perform similarly. Inside the bold ellipses, 
nonlinear scales are included, up to the innermost ellipse that corresponds to k mnx = 0.5 Mpc -1 . 
Here, f| n (i + 5) greatly outperforms Pg. See Ref. |3 1 for more details. 
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